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ABSTRACT 

We explore the th ermal state of the neu tron star in the Cassiopeia A supernova remnant using 
the recent result o flHo&Heinkel d2009l) that the thermal radiation of this star is well-described 
by a carbon atmosphere model and the emission comes from the entire stellar surface. Starting 
from neutron star cooling theory, we formulate a robust method to extract neutrino cooling 
rates of thermally relaxed stars at the neutrino cooling stage from observations of thermal 
surface radiation. We show how to compare these rates with the rates of standard candles - 
stars with non-superfluid nucleon cores cooling slowly via the modified Urea process. We 
find that the internal temperature of standard candles is a well-defined function of the stellar 
compactness parameter x = r g /R, irrespective of the equation of state of neutron star matter (R 
and r g are circumferential and gravitational radii, respectively). We demonstrate that the data 
on the Cassiopeia A neutron star can be explained in terms of three parameters: ft, the neutrino 
cooling efficiency with respect to the standard candle; the compactness x; and the amount of 
light elements in the heat blanketing envelope. For an ordinary (iron) heat blanketing envelope 
or a low-mass (< 1CT 13 M G ) carbon envelope, we find the efficiency ft ~\ (standard cooling) 
for x < 0.5 and ft ~ 0.02 (slower cooling) for a maximum compactness x « 0.7. A heat 
blanket containing the maximum mass (~ 10~ 8 M ) of light elements increases ft by a factor 
of 50. We also examine the (unlikely) possibility that the star is still thermally non-relaxed. 

Key words: dense matter - equation of state - neutrinos - stars: neutron - supernovae: 
individual (Cassiopeia A) - X-rays: stars 



1 INTRODUCTION 

Neutron stars are born with very high internal temperatures (T > 
10" K) but gradually cool down. Observing the thermal radiation 
from cooling isolated neutron stars and comparing their tempera- 
tures to theoretical models allows one to explore the (still uncer- 
tain) phys ical properties of superdense matter in neutron star cores 
(see, e.g.. IPethicMl 19921; iTsumtdl 19981: lYakovlev & Pethickll2004 
IPage. Geppert & WebeJl2006l : |Page et alj|2009l . for review). 

The supernova remnant Cassiopeia A (C as A), at a dis- 
tance of d = 3.4+°^ kpc teeed et all 1 19951) . is one of the 
youngest-known in the Milky Way. The supernova that gave 
rise t o the remnant may have been observed in 1680 dAshwortr] 
19801). though this may have been a mis-identification (see, e.g., 
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I2002L for review and discussion); nevertheless, a supernova ex- 
plosion in 1681 ± 19 years h as been estimated by examining 
the expansion of the remnant dFesen et al.l 120061) . Hereafter we 
shall assume that Cas A has a current age of 330 ± 20 yr. Al- 
though the supernova remnant is extremely well-studied, the cen- 
tral compact source was only identifi ed in first-light Chandra 
X- ray observations dTananbaumlll999h and subs e quently studied 
bvlPavlov et al.l J2000l) : IChakrabartv et alj ll200lh: IPavlov & Lunal 
d2009h . Recently, it was shown dHo & Heinkell2009h that the com- 
pact source in the remnant is a neutron star, which we will refer 
to as CANS (for Cas A Neutron Star or CArbon Neutron Star), 
with a carbon atmosphere and low magnetic field (B < 10" G). 
From the spectral fitting of the Chandra observations, the neu- 
tron star (gravitational) mass M, (circumferential) radius R, and ef- 
fective (non-redshifted) surface temperature T s were found to be 
M x 1.5 - 2.4M Q , R * 8 - 18 km (at 90% co nfidence), and 
T s ~ 2 x 10 s K, respectively dHo & Heinke] |2009|) . Prior to this 
last work, the possibility to study the CANS thermal evolution was 
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hindered by the assumption that the measured temperature reflected 
a local hot spot, hence it c ould only be used t o set an upper limit 
on the surface temperature dPavlov et al ] l20fB . However, with the 
recent results showing that the X-ray emission arises from the en- 
tire neutron star surface, T s can be treated as the average surface 
temperature. Because the next youngest neutron stars, for which 
surface thermal emission has been detected, have ages exceeding 
a thousand years, the CANS serves as a valuable window into the 
early life of a cooling neutron star. 

In this paper we use the theory of neutron star cooling and ana- 
lyze the thermal state of the CANS. Because cooling theory gives a 
wealth of cooling scenarios involving many uncertain properties of 
superdense matter in neutron star interiors (such as composition, 
equation of state - EOS, superfluid properties of baryons, etc.), 
we do not test these scenarios one by one. Instead, we propose 
a robust method to analyze observations of neutron stars of age 
~ 10 2 — 10 5 yr, i.e., those which have thermally relaxed interiors and 
cool via neutrino emission from their cores. We describe the proce- 
dure to extract the neutrino cooling rate of neutron stars (the ratio 
of their neutrino luminosity to the heat capacity) from observations 
and to compare the measurements with the theoretical cooling rates 
of neutron stars with non- superfluid nucleon cores which cool via 
the modified Urea process of neutrino emission. These theoretical 
cooling rates can be treated as standard neutrino candles, and the 
comparison between observations and theory can give valuable in- 
formation on the physical properties of neutron star interiors. We 
apply this method to interpret observations of the CANS, evalu- 
ate the CANS cooling rate and outline various physical models of 
CANS internal structure that are compatible with the observations. 
We also explore the possibility that the CANS is still in a thermally 
non-relaxed state. 



2 DATA ANALYSIS 

We summarize here the observations and spectral fitting of 
the CANS u s ed in our s tudy; further d etails can be found in 
iHo & Heinkel d2009h and iHeinke & Hoi ( bold) . Sixteen sets of 
Chandra X-ray Ob servatory archival data are considered (see 
IHeinke & Hoi l20ld for a listing), all using the ACIS-S charge- 
coupled device which provides spatial and spectral information 
jGarmirdl2003h . A series of observations, totalling 1 megasecond 
(the majority of a ll extant data), was taken in 2004 to study the su- 
pernova remnant dHwang et alj2004) . Several medium-length (~50 
ks) observations have been taken at various times between 2000 to 
2009. A 70-kiloseco nd observation in 200 6 was designed to study 
the compact source dPavlov & Lunall2009|), using a subarray con- 
fi guration to reduce ins trument pileup dDavisll200ll) . As discussed 
in IHeinke & Hoi d201Ct> . most of the 2004 dataset is subtly con- 
taminated by the presence of bad pixels affecting the flux from 
CANS. Although Chandra response functions include the effect of 
the bad pixels, the calibration becomes rather different, inducing 
small (< 2%) variations in the computed temperature, depending 
on telescope orientation (as the position of the CANS with respect 
to the bad pixels varies). We choose to include all data and per- 
mit variations in the fitted temperatures, while forcing all other pa- 
rameters to be the same between different observations (except the 
pileup grade-migration parameter a, permitted to vary with differ- 
ent frame times). We justify including all data by the much stronger 
constraints we obtain on neutron star mass and radius when we in- 
clude all the data, and by the fact that the temperature variations 
between observations are small (thus any effects on, e.g., the in- 



ferred neutron star mass and radius, will be <k 10%). CIAO 4.1 and 
XSPEC 12.4.0 are used for the data reduction and analysis. 

We note t hat using seven Chan dra ACIS observations from the 
past ten years . IHeinke & Hoi d2010h showed that the temperature of 
the CANS has been declining. However, the relative change is only 
ss 4% during the period of observation. Since we are considering 
the CANS in the context of its long-term temperature evolution in 
the present work, we treat the observed temperature as a constant 
at the current epoch. 

To fit the spectrum of the CANS, we used neutron star at - 
mosphere models with various compositions l lHo & Heinkell2009h . 
These models are constructed assuming a plane-parallel atmo- 
sphere (since the atmosphere thickness ~ 10 cm is much smaller 
than the stellar radius). The atmosphere is in hydrostatic and ra- 
diative equilibrium at constant surface gravitational acceleration 
g = GM/(R 2 VT^T), where 
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M 
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is the stellar compactness parameter, i.e., the ratio of its 
Schwarzschild radius, r g = 2GM/c 2 , to R. The efficient sep- 
aration of light and heavy eleme nts results in atmo s pheres 
composed of a single elem ent dAlcock & Illarionovl Il980l: 
iBrown. Bildsten & Ch ang 2002); the opacities are obtained from 
tables computed by the Opacity ProjecQ. Fur ther detail s of th e 
atmosphere model construction are given in IHo & Lail d200ll) ; 
lHo& Heinkel d2009l) . Only a non-magnetic carbon atmosphere pro- 
vides both a good fit to the spectrum of the CANS and an emis- 
sion size consistent with theoretical pr edictions for neutron star 
rad ii dHo & Heinkel 120091). In addition, IChang & Bildsten! d2004 
and lChang. Bildsten & Arrasld2010l) find that the high temperatures 
present in young (< 10 3 y) neutron stars remove all surface hy- 
drogen and helium on timescales shorter than the age of the neu- 
tron star. The (minimum) column depth and total mass of carbon 
required to form an optically thick photosphere is a 50 g cirr 2 
and ~ 10 15 g, respectively, and the density at one optical depth 
is ~ 0.02 - 4 g cirr 3 . Thus we only consider hereafter the carbon 
atmosphere model. 

Fits of the CANS spectrum by a neutron star atmosphere 
model depend on several parameters (in addition to corrections 
from pile-up and dust scattering). They are neutron star radius R 
and mass M, effective surface temperature T s , and effective column 
density of interstellar hydrogen Nh - In our spectral fits we have 
fixed the distance to d = 3.4 kpc dReed et al.|[l995l) . The small- 
est uncertainty is in N H (to within ~ 15% around /V H = 1.8 x 
10 22 cirr 2 ). The measured flux implies a bolometric (1 eV-10 keV) 
luminosity of 7.5+, \ x 10 33 erg s _1 at 90% confidence, assuming a 
distance between 3.3-3.7 kpc. 

Spectral fits give us allowed regions in T s , M and R param- 
eter space at different confidence levels. These results are shown 
in Figs. [T] and [2] The best-fit values are log T s [K] = 6.326, 
M = 1.65 M Q and R = 10.3 km (shown by filled squares). Fig. 
[T] presents 68% (dotted lines), 90% (short-dashed lines) and 99% 
(long-dashed lines) contours on the mass-radius diagram for neu- 
tron stars. The shaded upper-left corner is forbidden by the causal- 



ity principle. T he M — R diagram and causality constr aints are dis 
cussed, e.g., in lHaensel. Potekhin & Yakovlevl d2007l) . 
uncertainties are seen to be wide. 

To compare with theory, the solid lines in Fig.[T]give theoret- 
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Figure 1. (Color on line) CANS mass versus radius as inferred from spec- 
tral fits at 68% (dotted line), 90% (short-dashed line) and 99% (long-dashed 
line) confidence. The filled square gives the best-fit values. The shaded re- 
gion in the upper l eft corner is forbidden by the causality principle (e.g., 
lHaensel et alj|2007r) . The lower shaded region approximately shows the al- 
lowed values of M and R for realistic neutron star models. The three solid 
lines show M-R relations for neutron stars with the SLy, APR and PAL 1240 
EOSs. Filled dots indicate the most massive stable neutron stars for these 
EOSs. Triangles are the stellar masses above which the direct Urea process 
becomes allowed for the APR and PAL1240 EOSs (this process is impossi- 
ble in a stable star for the SLy EOS). 



ical M(R) relations for neutron stars which contain nucleon core s 

' iDouchin & HaensejJioOlh ; 



with three diff erent EOSs: SLy EOS by 



APR EOS by lAkmal. Pandharipande & Ravenhalll Jl998h: a ver- 
sion P AL 1 240 of the PAL EOS bv lPrakash. Ainsworth & Lattimen 1 
Jl988h . B y APR we mean the paramet r ization of APR results sug- 
gested bv lHeiselberg & Hiorth-Jensenl ll999|): sp ecifically, we use 
version APR I proposed bv lGusakov et alj d2005l) . The PAL 1240 is 
a version of the PAL EOS with the compression modulus of sat- 
urated nuclear matter K = 240 MeV and model 1 for symmetry 
energy jPrakash et alJI 19881) . Filled dots on the solid lines indicate 
maximum-mass neutron star configurations (M = M raax ). Triangles 
correspond to the mass thresholds M D u where the powerful direct 
Urea process of neutrino emission becomes allowed (more massive 
stars cool rapidly if the direct Urea process in their cores is not sup- 
pressed by superfluidity; see Sec. |4.2| (. For the SLy EOS we have 
A^max = 2.05 M Q , and rapid cooling is forbidden at any M < M max . 
For the APR EOS, M DU = 1.829 M Q and M raax = 1.929 M , while 
for the PAL1240 EOS, M DU = 1.25M Q and M max = 1.95M Q . The 
central shaded region in Fig.Q]shows the range of M and R that can 
be obtained for other EOSs which we consider as realistic. We see 
that the CANS limits on M and R, inferred from observations, are 
wider than the range of theoretical predictions. 

In Fig. [2] we present contours of 7* s inferred from the CANS 
observations as a function of R (left), M (middle) and x (right). 
These contours are given by the upper curves. They are projections 
of the appropriate T s - M - R regions (at 90% and 99% confidence 



levelsj^] on the corresponding axis, R, M or x. Having T s - M - R 
regions, we also show the associated - M - R regions, where 
j,oo _ yi _ x is the (redshifted) surface temperature as mea- 
sured by distant observer. The lower curves in Fig. [2] are the corre- 
sponding contours of Tf. Other curves are explained and discussed 
in Sectionl4~3l 



3 THREE COOLING STAGES 

Generally, neutron stars cool via neutrino emission from their inte- 
riors (mainly from the core) and via heat conduction to the surface 
and successive thermal surface emission. Cooling theory solves 
the heat conducti on problem within the star in General Relativ- 
ity jThorneiri977|) . accounting for internal neutrino energy losses 
and surface emission of photons. One can distinguish three main 
cooling stages. During the first stage (which lasts from ~ 10 
yr to a few centuries depending on the specific model), a newly 
born neutron star is thermally non-relaxed, with the core being 
cooler than the crust because of stronger neutrino emission in 
the core (seelNomoto & Tsurutalll98lUl987l;lLattimer et alj[l99l 
iGnedin. Yakovlev & Potekhinl 1200 it Ishternin & Yakovlevl 120081 
for details). The core is thermally decoupled from the crust; the sur- 
face temperature T s reflects the physics of the crust. During the sec- 
ond cooling stage (which lasts for 10 5 - 10 6 yr), the star is thermally 
relaxed inside; it cools mainly via neutrino emission. The main 
temperature gradient is located in a thin layer near the surface (in 
the so c alled heat-blanketing envelope, at densit ies p < pt, ~ 10'° 
g cm -3 . iGudmun dsson. Pet hick & Epsteinll 19831) . During the third 
(final) stage, the thermally relaxed star cools via the surface emis- 
sion of thermal photons. During each cooling stage, the thermal 
radiation of neutron stars carries different information on the stel- 
lar structure. Isolated neutron stars, whose thermal surface radia- 
tion has been detected, are mostly at the (second) neutrino cool- 
ing stage (with isothermal interior). The CANS is the youngest of 
these neutron stars and is thus of special interest. We study the ther- 
mal s tate of this neutron star using our general relativistic cooling 
code dGnedin et ai]|200ll) and useful semi-analytic expressions de- 
scribed below. 



4 NEUTRINO COOLING OF THERMALLY RELAXED 
STARS 

4.1 Basic relations 

A thermally relaxed star has an isothermal interior which extends 
from the center to the heat blanketing envelope (i.e., at p > pb). 
Takin g into account the effects of General Relativity (e.g.. lThornel 
1 19771) . isothermality means spatially constant redshifted internal 
temperature T(t) = T(r, t) exp(0(r)), where T(r, t) is the local inter- 
nal temperature, r is radial coordinate and O(r) is the metric func- 
tion that determines gravitational redshift. The cooling equations 
in such a star are much simpler than the general equations of heat 
transport. They reduce to the equation of global thermal balance 
dGlen & Sutherland! 1980h . 



2 Note that contours in Fig. \T\ refer to two-dimensional (M - R) confi- 
dence levels, while contours in Figs.[2J2]are for somewhat different three- 
dimensional (T s — M — R) confidence levels. However, this difference does 
not affect our principal conclusions. 
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Figure 2. (Color on line) Non-redshifted and redshifted effective surface temperatures T s and T™ (upper and lower lines, respectively) of the CANS, inferred 
from observations at 90% (short-dashed lines) and 99% (long-dashed lines) confidence levels, as a function of radius (left), mass (middle) and compactness 
(right). The shaded regions are theoretical values of T s and r™ for neutron star standard candles (having non-superfluid nucleon cores and cooling via the 
modified Urea process) with realistic M and R from the shaded region in Fig.[T] The solid lines refer to neutron star models constructed using the SLy, APR 
and PAL1240 EOSs with the direct Urea process switched off. All theoretical surface temperatures are calculated assuming iron heat blanketing envelopes, 
except for the two lightly-shaded regions in the right panel which assume carbon envelopes of mass AM = 1CT 8 M Q . 



C(T) ^ = -L?(T) - L?(T S ), 



(2) 



where L™(T) and L"(T S ) are, respectively, the neutrino luminosity 
and the photon thermal luminosity of the star (redshifted to a distant 
observer), while C(T) is the stellar heat capacity. These quantities 
are given by 



dVQ v (7»exp(20(r)), 



L7(T) = 
C(f) = j dVC v (T,p) 



AnaT^R 2 (1 - x) , 



(3) 
(4) 
(5) 



where Q v {T,p) is the neutrino emissivity, Cv(T,p) is the spe- 
cific heat capacity, cr is the Stefan-Boltzmann constant, and dV = 
4nr 2 dr exp(/i) is the element of proper volume determined by the 
appropriate metric function A(r). 

The hydrostatic stellar structure [p(r), O(r), etc.] is calcu- 
lated by solving the Tolma n-Oppenheimer-Volkoff equation (e.g., 
IShapiro & Teukolsky|[l983h . Then the cooling problem can be eas- 
ily solved by integrating l|2}; the solution gives T(t) and T s (t). 
To this aim, one needs to calculate the two functions, L™(T) 
and C(T) [by integrating over the stellar volume in $3$ and 
l[5}]. In addition, one needs to relate T s to the temperature Tb 
at the bottom of the heat blanketing envelope (p = p b ) by 
solving separately the stationary one- dimensional heat conduc- 
tion problem within this envelope (e.g . .lGudmundsson et aljfl983l ; 
IPotekhin. Chabrier & Yakovlev! 1 1997|- IPotekhin et alj 120031) . Ap- 



proximately, T s oc T^' z dGudmundsson et all 19831) . Since the en- 
velope is thin, one can set exp(C>(r)) a exp(0(/\)) = Vl - X within 
it, so that T a T b Vl - x. Eq. (O is solved with some initial con- 
dition, T(0) = T . However, for realistic values T > 10 9 K, the 
dependence T(t) at t > 1 yr is insensitive to To (the initial condition 
is forgotten at such t). 

Thus the cooling of thermally relaxed neutron stars is gov- 
erned by the functions L"(T), C{T) and T s (T b ), which are the only 
three functions that can be tested by comparing neutron star cooling 
theory with observations. This cooling theory is fairly insensitive to 



the microphysics of the stellar crust (at p > p b ) and to the thermal 
conductivity in the isothermal interior. The functions L™(T) and 
C(T) are sensitive to the (uncertain) microphysics in the neutron 
star core, while the r s (7b)-relation is based on the (better-known) 
plasma physics in the outer neutron star envelope. 



Our cooling model implies that the isothermal stellar inte- 
rior (p p b ) is spherically symmetric. The structure of the heat 
blanketing envelope (p < p b ) can deviate from spherical symme- 
try under the effects of strong magnetic fields. If the latter effects 
are substantial, T s varies over the neutron star surface, which is 
taken into account in our envelope m odels dPotekhin etai]|2003l : 
IPotekhin. Chabrier & Yakovlevll2007t) . However, we have checked 
that the magnetic field B < 10" G, which can be present in 
the surface layers of the CANS, does not affect L™, averaged 
over the surface. Thus, we set B = 0, so that the envelope 
models are determined only by the surface gravity g and by the 
mass AM of light (accreted) elements. This assumption is sup- 
ported by the non-de t ection thus far of pulsations from the CANS 
jMurrav et all |2002|; iMereghetti . Tiengo & Israel 200^; iRansoml 
l2002UPavlov & Lu na 2009; Halpern & Gotthelf 2 O10h . 



Our basic model will be the standard iron heat blanketing en- 
velope. The presence of a small amount of carbon (AM ~ 10~ 18 M o ) 
in the CANS atmosphere does not affect the cooling properties 
of the star. We also consider models with more substantial light- 
element envelopes. Light elements increase the thermal conduc- 
tivity and make the envelope more heat transparent (increase T s 
for the same T b ). Since we assume a carbon atmosphere, the enve- 
lope can contain carbon or heavier elements (hydrogen or helium, 
if present, would flow up to the surface). To this aim, we have con- 
structed new (carbon-iron) models with different carbon mass AM, 
and iron under the carbon layer. The thickest carbon layer, with 
AM ~ 10~ 8 M o , can extend to the density of a few times of 10 9 
g cirr 3 ; at higher densities carbon is destroyed by beta captures 
and pycnonuclear reactions. 
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4.2 Neutrino cooling rate 

The cooling problem is further simplified during the neutrino cool- 
ing stage. In this case can be neglected in lO, so that 



dT - - L°?(T) 
— = -i(T), t{T) = ^-J-. 
dt C(T) 



(6) 



The internal temperature is controlled by the physics of the stel- 
lar interior (mainly, the core), being insensitive to the structure of 
the heat blanketing envelope. The cooling is governed by the one 
function f(T), which specifies the cooling rate - the ratio of the 
neutrino luminosity to the heat capacity [K yr~'] - that is mainly 
determined by the properties of superdense core. It is only l(T) 
that can (in principle) be extracted from measured values of T s (t). 
However, the extraction is not straightforward. Two stars with the 
same mass and internal structure have the same T(t) (look the same 
from inside). However, if their heat blanketing envelopes are differ- 
ent, they have different surface temperatures T s (t) (look different 
from outside). 

We discuss the extraction, neglecting the effects of very fast 
rotation, superstrong internal magnetic field and possible devia- 
tions from beta-equilibrium. Then there is a single unique cool- 
ing function £(T) for the real EOS and microphysics of neutron 
star cores. This unknown function depends only on T and M (or 
some other quantity, R or x, instead of M). Theoretical models 
give different functions C(T), for different compositions, EOSs, and 
superfluid properties of the core. Of the two ingredients of £(T) 
[L™(T) and C(T)], the neutrino luminosity L™(T) is much more 
sensitive to the microphysics of superdense matter, than the heat ca- 
pacity C(T). Therefore, t(T) mostly reflects the neutrino luminosity 
of neutron star cores. The standard neutrino emis s ion mechanism 
is the modified Urea process (e.g., |Pethicklll992l ;l Yakovlev et al.l 
l200ll : jPage et alj 120061 . l2009h in non-superfluid cores. It gives 
L" MU (r) ~ (10 30 - 10 31 )r 8 8 erg s" 1 (where T g = T/10 8 K). We 
consider this as a neutrino standard candle (standard slow neutrino 
emission). The heat capacity of such cores is mainly determined by 
neutrons, C(T) ~ 10 38 Tg erg K . The cooling rate of the standard 
candle (SC) can be estimated as 



£ sc (T)~(03-3)Tl Kyr 



(7) 



We discuss (sc(T) in more detail in Section |473l 

The actual (unknown) cooling rate can strongly differ from the 
standard one, either enhanced or reduced. In this respect it is useful 
to introduce the ratio 



ft = e(T)/e sc (T), 



(8) 



which expresses the actual rate £(T) in units of standard candles (of 
the same M and R). It measures the neutrino cooling efficiency, the 
very important property of superdense matter to be extracted from 
observations of cooling neutron stars. The parameter ft should be 
taken at a temperature T(t) of the cooling star. If C(T) and £sc(T) 
have the same temperature dependence, then ft is just a number, 
independent of T; otherwise it depends on T, as discussed next. 

The enhancement of £(T) over {$c(T) can occur if more ef- 
ficient neutrino emission mechanisms are allowed in the inner 
cores of massive neutron stars. These fast mechanisms lead to neu- 
trino luminosity L™(T) oc T 6 . The strongest enhancement over 
L™ MU (with ft ~ 10 6 - 10 7 at f ~ 10 s K) can be provided by 
the direct Urea processes in nucleon or nucleon/hyperon matter 
dLattimer et alj[l99ll ; IPrakash et alJI 19921) . They open at densities 
higher than some threshold density p DU . Thus the neutrino lumi- 
nosity in non-superfluid low-mass neutron stars (where the direct 



Urea process is forbidden) is our neutrino standard candle, while 
at larger M it is much higher, leading to fast neutrino cooling. The 
density threshold p DU and the associated neutron-star mass thresh- 
old Mdu for fast neutrino cooling are very model dependent (see 
Sec, Hand Fig.0. 

Even if direct Urea processes are forbidden in the inner neu- 
tron star core, the neutrino luminosity can still be enhanced over 
L™ MV and can scale as Lf(T) oc T 6 . This happens provided the 
inner core has an exotic composition (contains pion condensates, 
kaon condensates or free quarks) owing to specific direct Urca- 
like processes in exotic matter. T he neutrino luminosity L"(T) 
depends on this composition (e.g., IPethickl [l992l : lYakovlev et al.1 
1200 ll ; IPage etai]|2006l [20091) but is smaller than L" DU (f), lead- 
ing to a smaller enhancement factor f c . In pion condensed matter 
one gets ft ~ 10 3 - 10 6 , while in kaon-condensed or quark matter 
ft ~ 10 2 - 10 4 . Calculations of neutrino emissivities and t(T) for 
these cases are very model dependent. 

The situation is even more complicated in the presence of su- 
perfluidity of baryons (Cooper pairing due to attractive component 
of the baryon-baryon interaction) in neutron star cores. In the sim- 
plest case of nucleon matter, one deals with superfluidity of neu- 
trons and protons. Hyperons and/or quarks (if available) can also be 
in a superfluid state. Pion or kaon condensates affect superfluidity 
of baryons. Each superfluidity is characterized by its own critical 
temperature T c t hat depends on p. Calculatio ns of T c (p) are model 
dependent (e.g.. lLombardo & Schulzell200ll) : various models give 
a large scatter of T c (p) values (~ 10 s - 10 10 K or higher). In central 
regions of massive stars, superfluidity disappears [T c (p) — > 0] be- 
cause short-range nucleon-nucleon repulsion destroys Cooper pair- 
ing. 

Baryon superfluidity can greatly mod ify neutrino emission , 
and hence affect i(T) and f c (see, e.g., lYakovlev et al.1 1200 ll ; 
IPage et aljfeood 120091) . First, strong superfluidity (T c » T) ex- 
ponentially suppresses the traditional neutrino reactions, which in- 
volve superfluid particles, and reduces C(T) because the gap in the 
energy spectra of superfluid particles blocks the reactions. Second, 
superfluidity initiates a different neutrino emission process due to 
Cooper pairing of baryons, which increases t{T). 

Let us illustrate these statements by taking neutron stars with 
nucleon cores as an example. Consider a low-mass star, which 
would cool via the modified Urea process in a non-superfluid 
case. Strong (for instance, proton) superfluidity in the core sup- 
presses the modified Urea process and reduces lij). If neutrons 
were non-superfluid, the main neutrino emission would be provided 
by neutrino pair production in neutron-neutron collisions (neutron- 
neutron bremsstrahlung); L™(T) would be proportional to T g (as for 
the modified Urea process) but would be weaker than L™ MV (T); the 
star would cool slower than its non-superfluid counterpart, leading 
to ft ~ 0.01 - 0.1. Now assume a mild triplet-state neutron super- 
fluid in the core, with maximum T c (p) of a few times of 10 s K. 
When the temperature in the core drops below this maximum, 
Cooper pairing of neutrons can initiate a sufficiently strong neutrino 
emission, with L™(T) oc T 8 . Such a star would be a faster neutrino 
cooler than the standard candle (with maximum ft ~ 10 - 100). 
On the other hand, this extra neutrino emission is negligible if the 
direct Urea process is open (even slightly) in the core. 

The effect of superfluidity in a massive star can be dramatic. 
The direct Urea process can be formally allowed in its core but 
suppressed by superfluidity. In this case, superfluidity transforms 
the fast neutrino cooling into a slow one and effectively increases 
M DU . It can greatly smooth out the transition from slow to fast cool- 
ing with increasing M. However, very massive stars would cool 
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rapidly (if fast cooling is allowed by the EOS) because superfluid- 
ity disappears in their central regions (see above). 

Recall that £(T) depends also on the heat capacity C(T). In 
a non-superfluid star with a nucleon core, C(T) is mainly deter- 
mined by neutrons (see above). The proton contribution is ~l/4 
of the ne utron one, and the contribution of electrons and muons 
is -1/20 dPagdl 19931) . Strong neutron superfluidity suppresses the 
neutron contribution, and strong proton superfluidity suppresses the 
proton contribution, but the electron and muon heat capacities al- 
ways survive. 

All in all, current theories predict different cooling rates £(T) 
which depend on many uncertain properties of superdense matter. 
The temperature dependence £(T) is mainly expected to be a power- 
law, 



£(T) = qT" 



(9) 



with n = 6 for fast neutrino cooling and n = 8 for slow cooling; q 
determines the cooling efficiency. From ([9} and ® at the neutrino 
cooling stage (after the initial conditions are forgotten), one obtains 
the cooling relation 



T(t) = [(n-2)qt]-' ,i "- 2 \ 
with 

1(f) = T/[(n - 2)t]. 



(10) 



(11) 



Such solu tions have bee n presented in previous works many times 
(see, e.g JPefhickll 19921) . If one knows /, n and T, then jl lb imme- 
diately gives the cooling rate t(T). For a fixed t the rate is directly 
proportional to T. Eq. J X Oi l greatly simplifies the solution of the 
cooling problem for a thermally relaxed star at the neutrino cooling 
stage. Note the weak dependence of T on t and on the cooling effi- 
ciency q, which stems from the strong temperature dependence of 
the cooling rate. 

It is instructive to introduce the parameter 



fr 



T(t) 



n-2 



((T) 



T S c(t) "sc - 2 { SC (T SC ) 



(12) 



where n and T(t) refer to a given star, while «sc = 8 and 7sc(?) 
refer to the standard candle of the same age, M and R. fr mea- 
sures the ratio of internal temperatures of the star and the candle, 
as well as the ratio of their cooling rates taken at their own internal 
temperatures. Because of the strong T-dependence of the cooling 
rate I [see eq. (5)], the parameter f T [oc l/^sc(7sc)] is drastically 
different from the ratio f e [oc 1/£ SC (T); see eq. ^8)]. According to 
eqs. l[9} and d 1 01 >. the cooling efficiency f [calculated at a given 
stellar temperature T(t)] is related to f T by 



ft = —T" 

qsc 



"sc - 2 1 

n-2 fsc- 2 ' 



(13) 



For instance, if the star has n = 8 and is twice as hot as the corre- 
sponding standard candle should be (f T = 2), its cooling efficiency 
is much lower than that of the candle (f ( = l/2" sc ~ 2 = 1/64). We 
will show in Sec tion |4T4l that f T can be easily constrained from ob- 
servations. However, it is f t that determines the efficiency of neu- 
trino emission and is the quantity of primary importance. 

We will also show that observations of the CANS are consis- 
tent with ft ~ 0.01 - 100. In this case it is reasonable to choose 
n = 8 (as we will do throughout the paper); f f and f T are then in- 
dependent of t and T. If in actuality n = 6, while we have chosen 
n = 8, then eq. dl3b underestimates ft by a factor of 1.5. 



4.3 Standard candles 

We simulate cooling of many standard candles, taking 17 nu- 
cleonic EOSs that describe the neutron star core. The EOSs 
includ e all 9 original ve rsions of the phenomenological PAL 
EOS ( Prak ash"et al J [l988l) . as well as 3 other versions of this 
EOS with the symmetry energy of nucl ear matter propose d by 
IPage & Applegatd ill 9921); the APR E OS l lAkmal et al .111 9981) ; the 
SLv EOS l lDouchin & Haensell 1200 ll) ; and three new modifica- 
tions (i-iii) of the PAL EOS. These modifications correspond to 
(i) Kp = 400 MeV and model 2 for the symmetry energy from 
IPrakash et al] Jl988h ; (ii) K = 300 MeV and mod el 3; and (iii) 
K = 400 MeV and model 3; the parameter tr in IPrakash et al .1 
i 19881) was taken to be zero in all three cases. These three EOSs 
are definitely superstiff and unrealistic. In particular, they predict 
the existence of very massive stars (with maximum masses up to 
~ 2.8 M G ). We include them to check that the scaling relations pre- 
sented below can be extrapolated to supermassive neutron stars. For 
any chosen EOS, we consider neutron star models of different M, 
from 1 M Q to the maximum mass, with a step size of 0.1 M e or 
smaller. For those models where the direct Urea process is allowed, 
we switch the process off to ensure that we consider only standard 
candles. In this way, we take about 170 neutron star models which 
cover a large area of the mass-radius diagram, with R from about 
9 to 15 km. We simulate cooling of all our models and determine 
T(t) and T s (t). In the simulations, we assume non-magnetic iron 
heat blanketing envelopes. However, the values of T(i) are almost 
insensitive to the composition of the envelopes (which determine 
T s ) at the neutrino cooling stage (see Seciion l-TTl i. 

We find that the internal temperature T(t) = Tsc(f) of all ther- 
mally relaxed neutron stars (standard candles) at the neutrino cool- 
ing stage can be approximated by 



r sc (f) = 3.45 x 10 8 K (1 - x) 



1+ °- 12 (l0km) IT) ' (14) 



where the compactness x is given by 10 and t c is some fiducial 
(normalization) timescale. The latter is set equal to the age of the 
Cas A supernova remnant {tc = 330 yr) that is convenient for our 
purpose. The time dependence of r S c in d 14b is taken in accordance 
with eq. dlOb . Calculations show that this time dependence is, in- 
deed, accurate at the neutrino cooling stage (see Section |5j- It first 
establishes in the core which thermally equilibrates faster than the 
inner crust. After the global internal thermal relaxation is over, it 
becomes accurate in the entire star (except for the heat blanket- 
ing envelope). The numerical factors 3.45 and 0.12 in d 1 4b are ob- 
tained by fitting the calculated values T(t c ). The accuracy of the fit 
is sufficiently good. The root mean square relative fit error is sa 3% 
(which translates into the 1.5% error of T s because T s oc 7* b 1/2 ); the 
maximum error over all models is ss 10% for an unrealistic (super- 
massive) model with M = 2.85 M Q ; the maximum error over ~160 
models with M < 2.4 M Q is « 7% (for the 1.9M Q model with the 
SLy EOS). 

Note, that our cooling curves contain systematic uncertainties 
which come from uncertainties in the calculations of the neutrino 
luminosity L^ MV and heat capacity C of the core. The main uncer- 
tainties are those in the effective masses of neutrons and protons, 
w* and m*, and in the squared matrix element |A4| 2 of the modified 
Urea process in the core. These quantities should be calculated self- 
consistently for a specific EOS (using specific many-body theory) 
as a function of density p, but such calculations (which would be 
most desired!) are currently unavailable. Instead, we set somewhat 
arbitrarily m* = m* = 0.7/n n (m„ being the free neutron mass), 
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R (km) M/M Q Compactness 

Figure 3. (Color on line) Cooling rates I [defined by eq. (6)] inferred from the observations of CANS at 99% confidence level (dashed lines) as a function 
of radius (left), mass (middle) and compactness (right). The rates are inferred assuming either an iron heat-blanketing envelope (AM = 0) or envelopes with 
outer carbon shells of mass AM = 10~ 12 M Q and 10~ 8 M o . On the right panel we add the dashed lines for AM = 10~ 13 and 10~" M Q . The shaded regions 
are theoretical £sc values for standard candles with M and R from the shaded region in Fig. [T] Solid lines refer to standard candles using the SLy, APR and 
PAL1240 EOSs. The right vertical scale on the right panel shows the internal temperature of the CANS. 



and we use th e stand ard approximate expression for |A4| 2 from 
lYakovlev etarl < l200ll) . To estimate the above uncertainties, let us 
introduce a typical nucleon effective mass mj, in dense matter and 
note that L™ MV oc (mJ,) 4 |A4| 2 and C(T) oc mj,. Then from eq. ( TTOb we 
conclude that T oc l/q [/6 oc (m^) -1 ' 2 1 All -1 '' 3 . For instance, allow- 
ing a 20% uncertainty in and 30% (uncorrelated) uncertainty in 
|A4| 2 , we obtain ~ 15% uncertainty in T and 7% in T s . We will not 
take these uncertainties into account in our further analysis (which 
is a common practice in the neutron star cooling theory) but warn 
about them. Taking the same nucleon effective mass for neu- 
trons and protons would be inaccurate in advanced simulations but 
is sufficient for estimating the uncertainties of T. 

Eq. J 14b gives a very simple method to calculate the standard 
slow cooling of thermally relaxed neutron star of any reasonable 
mass and radius during the neutrino stage. It does not require de- 
tailed knowledge of the internal stellar structure and microphysics 
of superdense matter. Because T is independent of the properties of 
the heat blanketing envelope at the neutrino cooling stage, eq. d!4t 
is not biased by these properties but allows one to take them into 
account while calculating r s . One can easily determine T for given 
values of t, M and R, and find T b = T / Vl - x at the bottom of the 
heat blanketing envelope. Then it is a simple task to calculate T s 
using the T^T^ relation. 

For example, the shaded regions in Fig.[2]represent the effec- 
tive surface temperatures (T s and T - ™) of standard candles of age 
t = 330 yr. Their masses and radii are taken from the (realistic 
EOS) shaded region in Fig. Q] All of them assume an iron heat 
blanketing envelope; the exception is for the two lightly-shaded re- 
gions in the right panel - these are for heat blanketing envelopes 
containing AM = 10~ 8 M G of carbon. Note a striking feature of 
theoretical standard candle curves - they are described by a well- 
defined (nearly universal) function of compactness x, but not of M 
and R. This is a remarkable manifestation of General Relativity. 
The presence of the maximum amount of carbon AM ~ 1O -8 M 
in the heat blanketing envelope does not violate the existence of a 
well-defined dependence of T s and T™ on x, but raises the effective 
temperatures by ss 25% for standard candles of the CANS age. 

At x < 0.5 the redshifted surface temperature of all stan- 
dard candles with iron blanketing envelopes is nearly the same, 



x 1.4 MK. With increasing x, the redshifted temperature T™ 
decreases, down to a 1.1 MK for standard candles with the max- 
imum compactness x k 0.7; in contrast, the non-redshifted tem- 
perature T s increases. Thus compact standard candles would ap- 
pear colder for a distant observer but would be hotter for a local 
observer. The differing behavior of and T s is another mani- 
festation of General Relativity. The weak sensitivity of r™(f) to 
variations of M and R for standard candles was first noted by 
Page & Applegatd dl992I) and later discussed in the literature (e.g., 
Yakovlev & Pethickll2004l . and references therein); here we quan- 
tify this effect. 

The three solid lines in Fig. [2] show T s and T™ for standard 
candles using the SLy, APR and PAL 1240 EOSs (their M-R rela- 
tions are displayed in Fig. [I}. We plot these curves using eq. d 141 ) 
and assuming iron heat blanketing envelopes. These curves are sim- 
ilar to each other, not only as a function of x, but also as a function 
of M (because the three EOSs are alike). 

Theoretical cooling rates [eq. {9), along with eq. <14H for stan- 
dard candles of the CANS age are shown by the shaded regions in 
Fig. [3] The three panels give these rates versus R, M and x. The 
masses and radii of neutron stars are again taken from the shaded 
region in Fig. [T] and the solid lines are again the theoretical rates 
for standard candles with the SLy, APR and PAL1240 EOSs. All 
these rates are independent of the properties of the blanketing en- 
velopes. We see again that there is a well-defined relation between 
the cooling rate I and compactness x. 

Since £(T) unambiguously determines T (for n = 8), the right 
vertical scale in Fig.[3]can be used to calibrate the internal neutron 
star temperature T. The expected values of T are a few x 10 s K. The 
non-redshifted local internal temperature T is somewhat higher. 
The relation between T and T is similar to that between and T s 
in Fig. [2] The difference between T and T in the neutron star core 
is larger than in the crust because of stronger gravitational redshifts 
in the core. The core is hotter (T is higher) than the crust for a given 
isothermal interior (given T). 
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4.4 Determining the cooling rate 

We now propose a simple and robust method to determine the cool- 
ing rate £{T) in a thermally relaxed star at the neutrino cooling 
stage. Assuming some values of M and R, a model for the heat 
blanketing envelope, and taking a measured (constrained) value T s 
(or 7"™), we can find (constrain) T. With a given (or assumed) age 
f, the cooling rate I can then be determined (estimated) from dl It 
(with n = 8 as discussed above). 

The application of this method to the CANS is presented in 
Fig. [3] which shows values of t (left vertical scale) and T (right 
vertical scale) of the CANS versus R, M and x. These values are 
inferred from spectral fits of the CANS at the 99% confidence level 
(Fig. |2j. The accuracy of inferring T is limited by the accuracy of 
our models for heat blanketing envelopes; the accuracy of inferring 
i is also affected by the assumption that n = 8. The upper (dashed) 
lines of I and T are obtained using the neutron star models with 
iron heat blankets. Lower lines are for models with carbon shells 
of different mass AM. We again see that the values of C and T are 
well-defined functions of x (for a given AM). The inferred values 
of I and T of the CANS for realistic values of M, R, x and AM vary 
within the factor f T ~ 3. The advantage of Fig. [3] is that it gives 
absolute values of the cooling rate £, an important parameter of the 
CANS. 

We now take the inferred values of T, use eq. l[8j and calcu- 
late the cooling efficiency ft, which is another important param- 
eter of the CANS cooling. As discussed in Section |4~2l relatively 
small variations of T yield much larger variations of ft. This is 
demonstrated in Fig. [4] which plots ft as a function of R, M and 
x. Lower, middle and upper lines refer to the blanketing envelopes 
with AM/M Q = 0, lCT 12 and 1CT 8 , respectively. The figure directly 
compares the CANS neutrino cooling rate to that of the standard 
candle. 

Figures[2}|4]impose constraints on the CANS cooling rate. Re- 
call that we can specify the rate by the single parameter f (provided 
we fix n = 8). Our results suggest that f t (as well as t, T s , T™) is a 
nearly universal function of x and AM. If the CANS has low com- 
pactness x < 0.4 - 0.5 and possesses either an iron blanketing en- 
velope (AM = 0) or an envelope with a very low mass carbon shell 
(AM < 10~ 13 M ), we find f ~ 1 (i.e., standard candle cooling). 
For larger x but the same heat blanket, we obtain lower f t (very 
slow cooling). For the maximum compactness x ss 0.7, we find 
ft ~ 0.02 - 0.03. If the CANS has a massive carbon heat blanket, 
then a higher cooling efficiency is required. With the most massive 
carbon shell AM ~ 10~ 8 M G at x < 0.4 - 0.5, the CANS cools faster 
than the standard candle, with ft ~ 30 - 50. However, ft decreases 
with increasing x and reaches the standard level at x ss 0.7. Note the 
opposite AM-dependence of the cooling rate I (Fig. [3} and cooling 
efficiency ft (Fig.[4]l. The star with smaller f cools slower and has 
higher temperature (and hence higher £) than the standard candle 
of the same age [see eq. t!3H . 

We thus set robust restrictions on the cooling efficiency ft of 
the CANS as a thermally relaxed cooling neutron star. Unfortu- 
nately, the values of f e do not specify unambiguously the neutrino 
emission mechanisms within the star (Section [4. 2t and can be re- 
alized by different physical models for the interior of CANS. For 
instance, ft ~ 1 is consistent with the modified Urea process in a 
non-superfluid neutron star core. It is equally consistent with en- 
hanced neutrino emission that is partially suppressed by superflu- 
idity. Thus there is ambiguity in the theoretical interpretation. We 
can explain the CANS data assuming different x and AM, while 



the values of ft, inferred at fixed x and AM, can be consistent with 
many models of CANS interiors. 



5 CANS AMONG OTHER COOLING STARS 

Figure[5]plots the CANS r™ limits on the r s °° - t plane, as well as 
the data on other cooling isolated neutron stars, The CANS limits 
are those obtained from 99% confidence contours (Fig. [2j. Obser- 
vational data on othe r sources are taken from references cited in 
lYakovlev et all d2008l) ; iKaminker et all d2009h . The CANS is the 
youngest in the family of cooling neutron stars whose surface tem- 
peratures are measured (constrained) more or less reliably. 

The thick solid line in Fig. [5] demonstrates the cooling of a 
standard candle. For this case, we take a neutron star with M = 
1.4M Q , R = 12.14 km (x = 0.34), the APR EOS in the core and a 
heat-blanketing envelope composed of iron. The initial part of the 
cooling curve is almost flat, which indicates that the stellar interio r 
is not thermally relaxed dLattimer etal]|l994rGnedin et alfeOOll) . 
The rapid drop of the surface temperature signals the end of thermal 
relaxation. By the current age of the CANS (« 330 yr), the star has 
an isothermal interior, and the results of Section|4]apply. The curve 
is reasonably consistent with the CANS observations, in agreement 
with the results of Section |4T4l 

The thick long-dashed cooling curve is calculated using the 
approximate cooling relation given by eq. d 1 4b . It is valid during 
the neutrino cooling stage in a thermally relaxed star. It accurately 
describes the cooling at t > 250 yr (when the thermal relaxation 
is over) until t < 10 5 yr (while photon cooling is still unimpor- 
tant). The densely-shaded region around this cooling curve is also 
calculated from eq. d!4t . We only show it in the limited range of t 
because eq. d 14b is inapplicable for very young and very old stars. 
This region is covered by cooling curves of standard candles (with 
iron envelopes); masses and radii are taken from that part of the 
central shaded region in Fig.Q]where x < 0.5. These cooling curves 
nearly coincide and are almost consistent with the CANS data (in 
agreement with the results of Section l4~4l see Figs. [2] and [3). The 
lower lightly-shaded region in Fig. [5] is covered by similar cooling 
curves of more compact standard candles (0.5 < x < 0.7) from 
the shaded region in Fig.Q] With increasing compactness, standard 
candles have lower T™ and cannot explain the CANS data. Neu- 
trino emission would need to be reduced to be consistent with the 
observations. A reduction of the cooling rate to ft ~ 0.02 - 0.03 is 
required for the maximum compactness x x 0.7. 

The thin solid lines in Fig. [5] are again the cooling curves of 
the 1.4M G APR star (with the iron heat blanketing envelope). How- 
ever, in this case the efficiency of neutrino emission is artificially 
enhanced (the lower curve) or suppressed (the upper curve) by a 
factor of 30 (i.e., f ( = 30 or 1/30, respectively). The thin dashed 
lines are analytic approximations of these curves; they are calcu- 
lated from eq. d 14b . in which we introduce the factor fj in ac- 
cordance with eq. d 1 0| > to renormalize the efficiency of neutrino 
emission. There is a good agreement between the lower thin solid 
and dashed curves (when the star is thermally relaxed and cools 
via the neutrino emission). Also note that the 1.4M G APR star with 
ft = 30 cools almost as a very compact (x = 0.7) standard candle 
( ft = 1). On the other hand, we obtain a slight disagreement be- 
tween the upper solid and long-dashed curves (for our basic APR 
star with strongly reduced neutrino emission, ft = 1/30). This is 
because the cooling of neutron stars with very low neutrino emis- 
sion from the core can be affected by neutrino emission from the 
crust. When we artificially reduce the emission from the crust, the 
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Figure 4. (Color on line) Contours of CANS cooling efficiency ft [defined by eq. {8)] as a function of radius (left), mass (middle) and compactness (right) at 
the 99% confidence level. The heat-blanketing envelopes are composed of either iron (AM = 0) or carbon shells of mass AM = 10~ 12 M o or 10~ 8 M o . 
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Figure 5. (Color on line) Theoretical cooling curves 7* s °°(f) compared with 
observations of cooling isolated neutron stars. The thick solid line is the 
standard candle (the 1.4M Q APR star) with an iron heat blanketing enve- 
lope. Upper and lower thin solid lines are for the same star but the neu- 
trino emission rate is artificially reduced and enhanced, respectively, by 
a factor of 30. The thick long-dashed line is the same as the thick solid 
line but calculated from the analytic expression d I4> that is valid for ther- 
mally relaxed stars at the neutrino cooling stage. The densely shaded strip 
is covered by numerous cooling curves of standard candles having iron en- 
velopes and x < 0.5; lightly-shaded strip is the same but for standard can- 
dles with 0.5 < .V < 0.7. The upper and lower thin long-dashed lines are 
analytic approximations 1141 of the corresponding thin solid curves. The 
thin short-dashed line is the analytic approximation of the cooling curve for 
the 1.4M APR star with a blanketing envelope containing AM = 10~ 8 M o 
of carbon. CANS is the data point at ( x 330 yr, while the values for th e 
other data points are ta ken from references, cited in lYakovlev et al.l J2008h ; 
iKaminker et~ai] fe009l) . 



upper thin solid curve shifts to the upper thin long-dashed curve 
(the shifted curve is not shown in Fig. [5}- We have checked that 
this effect becomes noticeable only in neutron stars with very low 
ne utrino cooling rate (/; < 0.03). Simila r effects were discussed 
bv lKaminker. Haensel & Yakovlevl d200ll) who studied very slowly 
cooling neutron stars with the neutrino emission in the core greatly 
suppressed by strong nucleon superfluidity. Thus our analytic ap- 
proximation becomes inaccurate for ft <sc 0.03. One should also 
be careful in applying this approximation to the case of f ( » 30 
because it is reasonable to expect that the cooling rate index n can 
change from n = rise = 8 to n = 6 (see Section POt . 



Finally, the thin short-dashed line in Fig. [5] is the analytic 
cooling curve of the APR 1.4M Q standard candle with a (maximal 
AM - 10 _8 M o ) carbon heat blanketing envelope. Recall that during 
the neutrino cooling stage such a star has the same internal temper- 
ature as a star with an iron envelope; nevertheless, it has a hotter 
surface because the carbon envelope is more heat transparent. The 
thin short-dashed curve is close to the upper solid and long-dashed 
curves meaning that the surface temperature of a standard candle 
with the carbon envelope is nearly the same as that of a star with 
lower cooling rate (ft ~ 0.03) and iron envelope. If the CANS has 
a carbon envelope and is not too compact (x < 0.5), then its neu- 
trino cooling rate should be over that of the standard candle by a 
factor of 30-50 (in agreement with the results of Section l4~4l . As 
the compactness increases, the surface temperature of the standard 
candle with a carbon envelope will tend to the CANS range; for the 
highest compactness and a carbon envelope, CANS would cool as 
a standard candle. 

An analysis of current theories of is ol ated neutron star s 
is given, e.g., by lYakovlev & Pethickl (l2004h ; IPage et all j200ot) ; 
lYakovlev etal] d2008l) ; |Page et al.l J2009T) . We note that the surface 
temperatures of nearly all stars lie between the upper and lower 
thin curves in Fig. [5] (the observed isolated stars do not require 
extremely slow or fast neutrino cooling). Our results indicate that 
the coldest of them (like the Vela pulsar), which have usually been 
treated as rapid neutrino coolers, can in fact be very compact stan- 
dard candles. Note that the surface temperatures of many isolated 
neutron stars (especially the colder ones) can also depend on strong 
magnetic fie lds, which affects heat tran sport in the blanketing en- 
velopes (see lPotekhin et al.ll2003ll2007l) . 
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6 IS CANS RELAXED? 

As shown in Section [4] observations of the CANS are consistent 
with the assumption that the star is thermally relaxed and cools rel- 
atively slowly (with a cooling efficiency f t ~ 0.02 - 50). In this 
section, we examine the possibility that the CANS is cooling much 
faster but its surface is not too cold because it is not yet thermally 
relaxed. We stress that cooling theories with the standard micro- 
phys ics of neutron stars give relaxation time s shorter than 300 yr 
(e.g.. lLattimer et al ll994l ; lGnedin et al.l200ll ). This means that the 
possibility considered here requires non-standard cooling scenarios 
which delay relaxation. 

In a young neutron star, the crust is hotter than the core be- 
cause of lower neutrino emission in the crust. Relaxation implies 
thermal equilibration of the crus t with the core that is mainly regu- 
lated by the physics of the crust |Lattimer et alJl994l ; lGnedin et al .1 
l200lUShternin & Yakovlev! |2008[) . It can be delayed by lowering 
the thermal conductivity and neutrino luminosity in the crust and 
by considering neutron star models with a thicker crust (which in- 
creases the thermal diffusion timescale in the crust). 

Figure[6]gives a number of cooling curves T s (f) which are cal- 
culated for two neutron star models based on the PAL 1240 EOS. 
We assume that the neutron stars have non-superfluid cores (since 
we focus on delaying relaxation by tuning the physics of the crust). 
The first model has M = 1.86M Q and R = 12.9 km. Recall that the 
direct Urea process is allowed for this EOS at M DU = 1.25M .Here 
we do not switch this process off (as in the calculations presented 
in Section|4j. As a result, the neutrino cooling rate is much higher 
[with f( ~ 10 6 ; see eq. (8)] than for standard candles. The second 
neutron star model has M = 1.3M Q and R = 12.9 km. Its mass 
again exceeds Mdu, and the neutrino cooling rate is huge. This star 
has a thicker crust than the 1.86M Q star, and its relaxation time can 
be longer. The displayed CANS errorbars are estimates from Fig. [2] 

Curve 1 in Fig. [6] shows the cooling of the 1.86M Q star with 
the standard physics of the core and crust. The end of the ther- 
mal relaxation manifests itself in the drop of T s (t) at t ~ 40 yr. 
By the current age of CANS, the star is thermally relaxed; it is 
much colder than the CANS (in agreement with the results of Sec- 
tion |4j- The largest delay of relaxation can be achieved by low- 
ering the thermal conductivity in the crust. We de monstrate this 
by using a model thermal conductivity proposed by iBrownl fcOOOl) 
(curve 2 in Fig. [6j. This model corresponds to an amorphous crust 
and can be regarded as the lowest limit on the crustal conductiv- 
ity. The relaxation time is now longer than the CANS age. By tak- 
ing the thermal conductivity to be lower than the normal one but 
above the lower limit, we can match any point in the region be- 
tween curves 1 and 2 and thus explain the CANS data. However, the 
hypothesis of low conductivity contradicts observati ons of thermal 
relaxation in quasi-persistent X -ray transients (see IShternin et all 
|2007| ;|B rown & Cummin j|2009l , and references therein). 

Curve 3 shows the cooling of the 1.3M G neutron star (with a 
thicker crust) assuming standard physics in the crust. The relaxation 
time is longer than for the 1.86M G star (curve 1), but the delay is in- 
sufficient to explain the CANS observations. Curve 4 demonstrates 
the cooling of the same 1.3M s tar but with the lowest crustal ther- 
mal conductivity terownl 12000) , corresponding to curve 2. In this 
case, relaxation is rather slow. Again by taking intermediate val- 
ues of thermal conductivity, we can fit any point between curves 3 
and 4. 

Relaxation can be accelerated if the neutron heat capacity in 
the inner crust is suppressed by neutron superfluidity. This effect 
is illustrated by curve 5, which is calculated using the lowest ther- 
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Figure 6. (Color on line) Cooling curves calculated for 1.86M neutron 
star (solid lines) and 1.3M Q star (dashed lines) with the PAL1240 EOS 
compared with observations of the CANS. Curves 1 and 3 correspond to 
standard physics input. Curves 2 and 4 employ the lowest limit on the ther- 
mal conductivity in the crust. Curve 5 is the same as 4 but adding moderate 
neutron superfluidity in the crust. Curve 6 is the same as 5 but the star has 
most massive carbon shell in the heat blanketing envelope. 

mal conductivity but assuming moderate neutron superfluidity in 
the crust (c.f., curve 4). The relaxation time drops to 700 yr, about 
twice the age of CANS. Finally, curve 6 is the same as 5 but for the 
star with the carbon heat blanketing envelope (AM = 2 x 10~ 8 M o ). 
This star is the same inside but looks hotter from the outside and 
has the same relaxation time. 

Thus the observations of the CANS can be compatible with a 
neutron star that undergoes fast cooling (via the direct Urea pro- 
cess) but is still not thermally relaxed. However, this can only be 
done by employing a very low thermal conductivity throughout the 
neutron star crust, which seems unrealistic. 



7 CONCLUSIONS 

We have used the cooling theory of isolated neutron stars and sug- 
gested a robust method (Section [4.4t to infer the neutrino cooling 
rate I [defined by eq. l[6]l] of stars from observations of their ther- 
mal radiation at the neutrino cooling stage (ages t < 10 5 - 10 6 
yr) after the end of their internal thermal relaxation. We simulated 
cooling of many models (wide ranges of EOSs, masses and radii) of 
standard candles (neutron stars with non-superfluid nucleon cores 
which cool via the modified Urea process) and formulated simple 
relations (Section[43} which determine the cooling rates and cool- 
ing dynamics during the neutrino cooling stage. We have shown 
that ultra-compact standard candles with compactness x ~ 0.7 have 
noticeably lower redshifted surface temperatures r™ than their less 
compact counterparts (x < 0.5). Based on these results, we devel- 
oped a method to compare the inferred cooling rates to the rates of 
standard candles and to determine the neutrino cooling efficiency f( 
[given by eq. (8)] of observed neutron stars. We have shown (Sec- 
tion[4]l that many physical properties, which characterize the ther- 
mal state of a cooling neutron star (such as the surface and internal 
temperatures, cooling rate), are well-defined functions of the stel- 
lar compactness x, rather than the stellar mass M or radius R. Self- 
similarity properties in the cooling of standard candles [eq. J14H . 
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their universal dependence on the compactness x, and the relatively 
low r™ of ultra-compact standard candles have been previously 
overlooked in the cooling theory. The suggested method is practi- 
cal. Fundamental uncertainties in the physics of superdense matter 
in neutron stars are incorporated into the parameter ft (or the cool- 
ing rate €) which can be determined (constrained) from observa- 
tions. 

We applied our theoretical formalism to study the thermal 
state of the neutron star in the Cas A supernova remn ant. We 
base our study on the recent results of lHo & Heinkd j2009l) . which 
showed that the observed X-ray radiation from this neutron star 
can be interpreted as thermal radiation from a carbon atmosphere 
(and the radiation is emitted from the entire stellar surface). We 
discussed the data (Section [2j and provided a theoretical interpre- 
tation (Section[4j, assuming that the Cas A neutron star (CANS) is 
thermally relaxed. Our analysis shows that the CANS cooling rate 
is ( ~ 10 5 K yr 1 and its internal temperature T ~ (2 - 3) x 10 8 K 
(Fig- [3}. Its thermal state is mainly determined by the cooling ef- 
ficiency ft, compactness x and mass AM of light elements in the 
heat-blanketing envelope. The ma gnetic field B < 10" ' G expected 
in the surface layers of the CANS Ibotthelf & HalperrJ 2001 1 20091 ; 
iHo & Heinkell2009l : iFMpern & GotthelfeOld) does not affect its 
cooling. 

If the CANS has a moderate compactness x < 0.5 and pos- 
sesses an iron heat blanketing envelope or an envelope with a low 
mass shell of carbon (AM < 10" 13 M e ), we find ft ~ 1 (imply- 
ing standard candle cooling). If it has the same heat blanket but 
larger x we obtain a lower f t (slower neutrino cooling), with a min- 
imum ft ~ 0.02 - 0.03 for x raax a 0.7. For the most massive car- 
bon shell (AM ~ 10~ 8 M o ) and moderate x, the CANS would have 
ft ~ 30 - 50 (implying a cooling rate enhanced with respect to the 
standard candle), but for the same shell and maximum x it would 
have the standard-candle cooling efficiency fi~ 1. In summary, we 
find 0.02 < f t < 50 for the CANS. Though these results do not 
give an unambiguous physical picture of the thermal structure of 
the CANS, they do impose firm and robust constraints. 

We also examined the possibility that the CANS has a high 
neutrino cooling rate but is still not thermally relaxed and, there- 
fore, is not very cold (Section . We found that this is possible 
(though not likely) provided the neutron star has a very low ther- 
mal conductivity in the crust. 

Note that we considered heat blanketing envelopes made of 
iron, as well as envelopes composed of an outer carbon shell. We 
could have complicated the model, e.g., introducing an oxygen 
shell beneath the carbon shell, but we do not expect qualitatively 
different results. 

It would be instructive to apply the above analysis to other 
cooling neutron stars, in an attempt to reconstruct the cooling rate 
function 6. It is expected to be one and the same function for all 
neutron stars, and it contains an important information on neutron 
star structure. 
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